Ricci-Tersenghi et al. reply: In their comment Jl| to our paper |2| de Candia and Coniglio 
show evidence that the equilibrium overlap distribution P(q) of the frustrated Ising lattice 
gas (FILG), do not coincide with the one that could be estimated from the off-equilibrium 
correlation and response functions presented in |2[], for low temperatures and large chemical 
potentials. From the theoretical point of view one expects the relation between static and 
dynamic quantities to hold as long as one-time observables (like energy or density) converge 
to the equilibrium values [[§. 

In order to clarify what is happening in the glassy phase of the FILG, we have performed 
a new large set of simulations, working mainly at (3 J = oo and f3fi = 10. At equilibrium we 
find a P{q) with two clear peaks and a continuous part in between (like in 0). Whether the 
continuous part disappears in the thermodynamical limit, thus suggesting a one step replica 
symmetry broken (RSB) solution, is hard to decide at present. 

In the out of equilibrium regime (quenching from (3fi = — oo to = 10) we have repeated 
all the experiments presented in [Q], but now letting the density grow even after the waiting 
time t w . The results (see Fig. [I]) are indeed different from those where the density is kept 
fixed after time t w . In Fig. [I] we clearly see that for large waiting times the agreement 
between static and dynamic susceptibilities is very good indicating that the system is not 
trapped in any long-lived metastable state. 

Still the open question is why in the experiments presented in || (and confirmed in [|TJ) 
the responses are so different and seem not to be related to the corresponding static suscep- 
tibilities. The answer is related to the sudden change of dynamics at time t w : up to t w we 
evolved the system in a grancanonical ensemble where the density tends to increase thanks 
to the large chemical potential, while after time t w we kept the density fixed, allowing the 
system evolve only in a canonical ensemble. In order to verify this hypothesis we have done 
the following numerical experiment: after quenching the system we let it evolve for t w time 
steps within the grancanonical ensemble and for at w more in the canonical one. Finally 
at time (1 + a)t w we switch the field on and we measure the response [always with a fixed 
density p(t w )\. For a = we recover the behavior reported in [Q,^f. Nevertheless for large 



a the system "forgets" the drastic change of the dynamics done at time t w , thermalizing 
in the canonical ensemble at a fixed density p(t w ) and the dynamical susceptibilities are in 
agreement with the static ones measured at the same density as can be seen in Fig. |2|. 

It must be emphasized that, in both figures, the off-equilibrium results obtained with 
L = 40 systems are still compatible with linear fits and that the curvature of the static 
result (L = 10) may show strong finite size effects. 

Concluding, we have shown that the link between static and dynamic susceptibilities 
is valid on the time scales we can reach in our simulations for both the canonical and the 
grancanonical ensembles. Regarding the kind of RSB solution present in the glassy phase 
of the 3d FILG (whether it is characterized by one or infinite RSB steps), we believe that 
further studies are needed in order to clarify this point. The model still seems to be a good 
candidate for a 3d system with 1-RSB. 
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FIG. 1. Grancanonical response versus correlation. 
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FIG. 2. Canonical response versus correlation. 
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